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ABSTRACT 

There is a need for a slmplo theoretical approach to the calculation of Bun- 
light penetration depths suitable for passive multiapectral remote sensing of 
water resources. An earlier paper presented an approach which is readily 
adapted to this calculation and which provides reasonably good agreement with 
more accurate but time-consuming radiative transfer models. The needed 
modifications are described and the model is used to calculate the penetration 
of sunlight into clear ocean water at several wavelengths throughout the visible 
portion of the spectrum, calculations for both clear and turbid water are car- 
ried out for the two visible channels of the multispectral scanner on NASA's 
ERTS-1 satellite. The effect of a reflective bottom on the upwelling light field 
is discussed. Measurement parameters needed for the passive remote deter- 
mination of water depths are identified and the use of submerged reflective 
panels for surface truth measurements is discussed. 


2 


INTRODUCTION 

A number of Investigators are using the Multlspectral Scanner Subsystem (MSS) 
on NASA's Earth Resources Technology Satellite (ERTS-1) for remote sensing of 
water resources. 1 In many of these investigations It Is necessary to have a rough 
idea of the effective penetration depth of Imagery In each of the spectral bands 
employed by the MSS. Klemns has devised a slanting board with painted stripes 
chosen to match the two visible channels of the MSS and calibrated to give a 
visual estimate of this penetration depth/ Other Investigators use a white circular 
disk, called a Secchl disk, lowered into the water until It disappears from view In 
order to obtain approximate information as to the penetration of sunlight into 
the water. The Secchi disk has been widely used in oceanography for meas- 
urements of water transparency. Some ERTS investigators are attempting to 
use the upwelling light reflected from a shallow bottom to estimate water depths. 

In each of the above investigations, there is a need for simple theoretical 
techniques relating the observations to the optical properties of both the water 
and the submerged reflecting surfaces. In a recent paper 3 an optical model 
developed by Hansen 4 and Gordon 5 was used to calculate the upwelling spectral 
radiance emerging from the sea. Although the accuracy of the model is limited 
by the simplifying assumptions upon which it is based, it is superior to single- 
scattering models for calculations of ocean color spectra and is easily adpated 
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to the calculations needed boro. The purpose of this paper Is to Illustrate tho use 
of the moclol In estimating depths of sunlight ponotratlon In typical remote 
sensing situations. 

PENETRATION DEPTH CALCULATIONS 

With a few modifications of the equations presented In reference 3 we can 
obtain analytical expressions for the reflectance of the sea and tho apparent 
brightness of submerged reflecting panels as functions of depth. Skylight con- 
tributes less than 30% of the total upwelling radiance emerging from beneath 
the sea. We may, therefore, further simplify our calculations by neglecting 
this contribution and consider only the penetration of sunlight into the sea. The 
upwelling radiance predicted by the single-scattering model for this oase is 
given by equation (8) in reference 3. For a single, homogeneous layer of depth 
z in the sea, Illuminated by the sun at an angle of 9 from the nadir (the solar 
zenith angle is tt - 9), this equation may be written as 

_ N(M _ T 0 T(e) n (1 - «’>] (1 , 

1 H / V ~ 7f — 7TT - — ■ - — 

* 0 7J 2 (1 - COS 9 ') c(\) 

Equation (1) gives the nadir-viewed upwelling spectral radiance N(\ ) 
emerging from the sea divided by the incident solar irradiance H Q (\) at wave- 
length k. This ratio will be called the radiance reflectance, R(X), throughout this 
paper. T 0 is the transmittance of the air-sea interface for upwelling light, T is 
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the transmittance for downwolllng sunlight, II 0 Is the spectral irrndlancc on a 
horizontal piano due to Incident sunlight, 17 Is the Index of refraction of seawater, 

/* Is the spectral volume scattering function, and c Is the single scattering 
extinction coefficient. The quasi-single scattering extinction coefficient e* Is 
defined to be c(l - F), whore Is the single scattering albedo and F Is the 

forward scattering coefficient. Replacing c by c* in equation (1) converts It to 
tho quasl-single scattering case described in reference 3. O' is the angle in 
water made by tho rays of the sun with respect to the zenith. It is related to & 
by Snell's law: v sin {n - O') ~ sln(77 - r*). We are here assuming a perfectly 
smooth sea surface. 

We shall examine the predictions of eq. 1 for clear ocean water. In 
evaluating equation 1 we use the usual Fresnel reflectance formula for calculating 
the transmittances T Q and T. For 0 = 150°, as is approximately the case with 
ERTS-1 when it crosses the central latitudes of the U.S., we have T Q - .9796 and 
T(150°) = .9785. All of these quantities are defined in reference 3. 

The data wo shall use for c(\) and w 0 (\) in our simulation are taken from 
a paper by Tyler, Smith, and Wilson 6 for clear water optically similar to that 
found in the Sargasso Sea. Data for (3 is taken from measurements of the 
scattering function in the Sargasso Sea made by Kullenberg. 7 The index of 
refraction of the sea water, t) , is approximately 4/3 and is not strongly wave- 
length dependent in the visible portion of the spectrum. The forward scattering 
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coefficient F is calculated from by straightforward numerical integration. We 
shall detormlno the radiance reflectance using equation 1 at 450 nm, which is 
close to the wavelength of maximum transparency of the sea, and for a zenith 
sun. 

The results of this calculation, using Eq. (1), aro shown in Figure 1 for 
both the single-scattering (SS) case and for the quasi-single scattering (QSS) 
case (when c* is used for c in Eq. 1). 

We see from this figure that the light p jnotrates deeper and produces a 
stronger return signal with the QSS model than with the SS model. The depth for 
the radiance reflectance to x’each 90% of its maximum value is about 25 meters, 
according to the QSS model. We can also see from this figure the importance of 
including the effects of multiple scattering in calculations of penetration depths 
when is large. 

In order to extend these results to other wavelengths in the spectrum, the 
Kullenberg scattering data at 460 and 655 nm was linearly interpolated to give 
values for p at 450, 500, 550, 600, and 650 nm. In addition, the function at 350 
and 400 nm was set equal to that at 450 nm and at 700 nm it was set equal to that 
at 650 nm. The results for the quasi-single scattering model are shown in 
Figure 2. We see that the upwelling red light reaches a maximum quickly, in 
less than one meter of depth. The reflectance is correspondingly low for this 
color. Blue light penetrates much more deeply and returns a stronger signal. 
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From Figuros 1 and 2 wo can see the decreasing Importance of deeper and 
deeper layers to the upwelllng radiance spectrum. Stnco in reality the optical 
properties of the ocean vary with depth, when one measures an ocean color 
spoetrum he obtains a weighted average of those properties over depth, if the 
depth dependence of these properties is known, then It may be accounted for by 
replacing equation (1) with 


R 


i 

T T( ') f* r* 

(K) - .-(a , r 1 , z)exp[“ (1 " see t ') c* (z") dz"]dz'. 

Jj 2 COS >■' d 0 J 


( 2 ) 


To show more explicitly the wavelength variation of the penetration, the 
radiance reflectance calculations wore used to determine the layer depth for 
which the upwelllng radiance is 90% of its maximum value. The resulting so- 
called "penetration depths” are listed in Table 1. This data applies to clear 
ocean water only. Also given are the SS and QSS attentuatlon lengths l/c and 
l/c*, respectively, and the approximate horizontal visibility distance using 
Duntley's criterion. These results are plotted in Figure 3. 

ERTS PENETRATION DEPTHS 

We shall apply the above results to the ,/o visible- wavelength bands of the 
ERTS-1 MSS. The "green" band, covering the spectral range 500-G00 nm, 
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is called band 4 in the ERTS Data User's Handbook. The red band, 
covering 600-700 nm, is called band 5. In performing this calculation, 
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wo again used tho Tylor, Smith, and Wilson dntn 0 to calculate c(M and ^’ 0 (^) at 
20 nm Intervals over the range from 500 to 700 nanomotors. The scattering 
functions used In tho previous section wore linearly Interpolated to obtain i at 
tho same wavelengths used for c{k) and (X). The resulting radiance reflec- 
tances wore then averaged over tho two EllTS-1 MSS visible wavelength bands 
separately. Tho rosultB nro plotted in Figure 4 for a 30* zenith sun angle and 
for nadir viewing. 

We see in Figuro 4 that the SS model is adequate for bond 5 but the QSS 
model should be used in calculating band 4 penetration depths. The depths for 
90% of maximum radiance reflectance are approximately 18 and 3.8 meters for 
bands 4 and 5, respectively, using tho predictions of the QSS model. The low 
penetration of band 5 is due primarily to a high extinction coefficient for sea 
water over the wavelength range of this band. 

Since most ERTS-1 oceanographic investigations ore concerned with coastal 
waters, which are generally more turbid than open ocean waters, the above re- 
sults will be of limited usefulness, other than to indicate the clear water limit of 
sunlight penetration into the sea. In order to provide an estimate of the results 
one might obtain in a more turbid area, we turn to the measurements of fi and c 
made by Petzold at 515 and 520 nm, respectively, in the relatively turbid waters 
of San Diego harbor. 10 Although significant variations in tiie shapes of the 
wavelegnth spectra of fi and c are known to occur when going from clear to 
turbid water, lack of suitable measurements of these variations compels us to 
assume that the changes are negligible. Since we are mostly concerned with 


obtaining an estimate of the range of penotrntlon depths which might bo found 
with I3RTS-1 Imagery, the crudeness of our assumption should not bo too 
onorous. 

The total scattering coefficient b Is defined by the relation 

t> f ( y ) d < 2 ' f ,) s i n . d » 

*o 

► 

where >- is the scattering angle. Petzold obtained a value- of 1.824 for b at 615 nm. 
For c ho measured 2.19 at 530 nm, giving an approximate -. 0 of 0.833. The 
values for b, c, and .» 0 at 500 nm used in the above clear water simulation were 
0.03234, 0.0846, and 0.3823, respectively. Tho turbid (San Diego Harbor) - to - 
clear (Sargasso Sea) - water ratios of these quantities are 56.4, 26.8, and 2.18, 
respectively. By multiplying tho clear water values for r and c by these ratios 
at each of tho wavelengths from 500 nm to 700 nm used in the previous simula- 
tion, we obtain a crude representation of tho requin data for turbid water. 
Repeating the calculations that led to figure 4, but using our simulated turbid 
water data, yields the results shown in Figure 5. Tire corresponding QSS 90% 
penetration depths for bands 4 and 5 are 1.8 and 0.16 meters, respectively. 

As expected, the depth of penetration is very much less for turbid water 
than for clear ocean water. Although this calculation is only approximate, it 
should give at least an order-of-magnitude estimate for turbid estuarine water. 


0 

By using bettor data for c and /*, tho nccurncy of the estlmnto can bo greatly 
Improved, 


REFLECTION FROM A SHALLOW BOTTOM 

If the lower surface of the ocean Is partially reflecting, as would be the 
case with a shallow bottom or a submerged reflecting panel, then additional 
terms must be added to equation 1 to account for tho additional flux reflected 
from the bottom. If this surface Is a lambertlan rofloctor, with Its diffuse 
Irradlanco reflectivity given by r (\ ) (defined to bo tho ratio of upwolllng to 
downwalling spectral L radiance Just above the surface), and If we Ignore the 
multiple reflections between the bottom and the water-air Interface, then tho 
added contribution to the upwelling radiance reflectance above the sea can be 
written in the quasi-singlq scattering approximation as 




VW ~c*<\)( 1 -sec 0 9 ) z 

— r 0 (\)e 

•nr; £ 




where z is the bottom depth. Adding (3) to equation (1) (with c replaced by c+) 
gives the radiance reflectance for an ocean having a reflective bottom of 
irradiance reflectivity r Q at depth z: 


R(\) = 


T o T(0) 


1 ‘ 


c*(l - cos #' ) 


7t 


(4) 
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The limitations of tills expression must bo pointed out. Firstly, QSS ealcula- 
tions of the downwclllng, forward-scattered (and unscatterod) light field aro not 
expected to be as accurate, due to the nature of the QSS assumptions (given in 
ref. 3), as were the calculations of tho upwelling, backward-scattered light field 
discussed In die previous section. In the presence of a reflecting bottom, this 
downwclllng light field will be re-directed upwards and will combine with path 
radiance due to backscattering to produce tho upwelling radiance spectrum 
emerging from the water. Thus, any inaccuracy in the predictions of the down- 
welling light field will be present in predictions of the bottom-reflected portion of 
the upwelling light field (given by cq» (3)), Thus, we cannot expect as good on 
agreement between tho QSS model and more rigorous radiative transfer calcu- 
lations when tiie omerging light rofloctcd from tho bottom is a substantial frac- 
tion of die total upwelling light field. 

These observations have been confirmed by Gordon. Mo lias found that tho 
QSS model predicts a radiance reflectance due to tho bottom only which is 20% 
less than that predicted by his Monte-Carlo radiative transfer calculation, for 
the case % = 0.2, a bottom albedo of unity, and a bottom depth of one optical 
thickness. 11, 12 

Another source of error can be identified which is not accounted for In 
the above calculations. When the bottom is bright and the water is clear and 
shallow, some of the light reflected from the bottom will be reflected down- 
ward by the air-water interface and will be again reflected by the bottom. 

This contribution, and the error it produces, will be substantial whenever 


the bottom depth is small. 
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Duo to those limitations, equation (4) cannot be expected to give accurate 
predictions of the up welling radiances In the presenco of a shallow bottom will) 
high reflectivity, VVe can, howovor, use equation (4) to discuss in a qualitative 
way the various functional dependencies between the upwclling light field and 
the various subsurface parameters Influencing it. 

For example, let us oxamlne the predictions of this expression for passive 
remote determinations of bottom depth. If one wishes to determine z from 
measurements of 11 (a), using equation (4), then r Q (a.) and the depth dependence 
of both p (h,t ,*') and c (a ) would have to be measured as well. This require- 
ment greatly restricts the usefulness of the technique. But if several wave- 
lengths ur t*-od, and a few simplifying assumptions aro mado, then it becomes 
theoretically possible to determine the depth with substantially less information 
being required. For example, suppose we limit our measurements to areas 
where the water is homogeneous and is sufficiently clear and shallow that the 
light from the bottom is much stronger than that scattered upward by the water. 
In this case the first term in equation (4) may be neglected and we are left with 
equation (3). Now let us suppose that measurements of lt(\) are made at two 
wavelengths, designated \ t and \ 2 . Writing equation (3) for each of these and 
then dividing the one for A , by the one for \ 2 , and solving for z, yields the 
expression: 


7 . ~ 


1 

[c*(^j) - c*(\ 2 )](l - sec 6 ) 


In 




R(\) " 
r 0^2^ 


( 5 ) 
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Using this expression, wo only need to know the difference in quasi-single- 
scattering extinction coefficients and tho ratio of r to 11 at the two wavelengths 
in order to determine tho depth. The slnglo scattering version of tills expres- 
sion has been used by Polcyn and Lyzenga 'o estimate bottom depths in the 
Bahama islands using ERTS imagery. 13 

Although tho accuracy of Equation (5) is limited, as was discussed before, 
it can be used to identify sevoral potential sources for considerable error 
which would bo present in calculations based on similar approaches. For 
example, an undetected change in the shape of the bottom reflectivity spectrum 
would be erroneously detected as a change in the water depth. An unexpected 
increase in water turbidity could be incorrectly interpreted as a decrease in 
bottom dopth. 

The difficulty of finding water of sufficient shallowness and clarity that the 
first term of equation (4) may be neglected makes this approach to accurate 
passive r^'uote determinations of bottom depth rather unattractive, even if a 
rigorous radiative transfer calculation were used in place of the present approxi- 
mate analysis. 

By dropping the assumption that the water is shallow and clear enough to 
neglect the first term of equation (4), and by assuming that measurements will 
be made of c(\), /i (A , 0 '), and r Q (\) at several wavelengths for some point in the 
field of view of the remote sensor, then the more generally valid equation (4) 
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can bo used. Ilowovor, the potential errors producod by unaccounted for changes 
in bottom reflectivity and water turbidity away from tho surface measurement 
location will persist, and, in tho author's view, will sevorly limit tho usefulness 
of the technique for routine measurements of bottom depth with the kind of regu- 
larity and accuracy which is desired. Any attempt to make passive remote 
measurements of water depths must recognizo the need for a number of optical 
measurements at the surface and the limitations of accuracy which will be 
imposed by the absence of any of those measurements. 

REFLECTION FROM SUBMERGED REFLECTING PANELS 

Let use now use equation (4) to determine the depth dependence of the up- 
welling spectral radiance integrated over bands 4 and 5 of tho ERTS-1 MSS for 
the case of submerged reflecting panels of known spectral reflectance r Q (\). 

We shall assume that the panel reflectance is constant over the wavelength 
ranges involved. Using the clear and turbid ocean cases described earlier, 
and evaluating equation 4 with assumed bottom albedos r of zero, 0.1, 1, 10, 

20, 40, and 60 percent, we obtain the results shown in Figures 6 and 7. 

We must again point out the limitations of these results. Although they should 
be more accurate than those obtained using the single-scattering approxi- 
mation alone, the results do not fully and accurately account for the effects of 
multiple scattering within the medium and they totally neglect the multiple re- 
flections taking place between the bottom and the surface. Nevertheless, Figs. 

6 and 7 can be used to illustrate the general behavior of the radiance reflectances 


with increasing depth, At very great doptlis, the reflecting panels effectively 
disappear and we are left with the more accurate caso of an infinitely deep 
ocean described in refei'once (3), 

Looking at the results shown in Figures 6 and 7, we see immediately that 
the use of a reflective panel having a i oflectance of 20% or greater can signifi- 
cantly extend the apparent depth of penctratton of sunlight into the sea. Visual 
estimates of this penetration depth using refloating panels or colored Socchi 
disks will depend upon the visual contrast of the submerged target, which will 
in general appear’ to decrease In size with increasing depth. The results show* 1 
in Figs. 6 and 7, however, are based upon the implied assumption that the hori- 
zontal dimensions of the submerged panel are substantially greater than the 
depths at which it is used. It is therefore difficult to relate the results of Figs. 

6 and 7 to the disappearance of a panel from view. Thus, great care must be 
taken when using Secchi disks or reflecting panels in attempts to estimate depths 
of sunlight penetration application to remote imagery of water resources. 

The measurement of Secchi depth is subject to a number of extraneous 
influences, most notably the disruption of the image due to surface waves, the 
extent of obscuration of the sun due to atmospheric ha ze, clouds, or overcast, 
and variations in the visual acuity of individual observers. As a result, the 
Secchi depth measurement is extremely subjective and inaccurate, and can be 
regarded as little more than a qualitative estimate of water clarity. 


The inclined submerged reflecting panel, properly used, can bo a powerful 
surface truth toolmique in remoto sensing situations. If its reflectance spectrum 
r Q (\) Is known accurately, and if its horizontal dimensions are much greater 
than the depths at which it is used, then, with a suitably designed radiometer 
it should be possible to obtain the full depth dependence of the radiance reflec- 
tance R(\) approximated by equation (4). For very small values of z, to a first 
approximation the second term dominates and it might be possible to use the 
measurements of R(\) in this range to determine the QSS extinction coefficient 
c*. 

For increasing values of z, the first term in eq. (4) will have a growing in- 
fluence, The value of c* obtained from the measurements at small values of z 
can be substituted into this equation and an estimate of the volume scattering 
function (i may be thereby obtained at the wavelength and scattering angle in- 
volved. The technique requires accurate knowledge of the down welling sunlight 
spectral irradiance and it cannot be used if the sky is hazy or overcast. The 
presence of skylight can significantly influence the results, especially at 
shorter wavelegths, and should be carefully considered in any intended appli- 
cation of the technique. 

CONCLUSIONS 

We have shown the utility of the quasi-single scattering model for calculating 
the penetration of sunlight into natural waters. Although the accuracy of the 
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model Is still not optimal, It is superior to single-scattering models for those 
calculations. The lack of accurate extinction and scattering data is at present 
the primary factor limiting the accuracy of the approach when applied to 
practical situations. 

For more accurate calculations of the upwelling radiance spectrum emerging 
from natural waters, or of the penetration depth spectrum, one must have the 
full wavelength dependence of both the extinction coefficient and the volume 
scattering function (over its full angular range) for the body of water being 
modeled, regardless of the particular theoretical model being used. This 
information can come from measurements made in the field or it can be of a 
suitably formulated microscopic model of the optical properties of the water, 
as described in ref. 3. It is essential that progress be made either in developing 
the mici’oscopic model or in providing the instrumentation needed to gather 
the information directly if accurate upwelling radiance spectra and penetration 
depths are to be calculated. This conclusion applies as strongly when a shallow 
bottom is present. And in this case it is also necessary to have the spectral 
reflectance of the bottom. 

When the water color spectrum is being modeled in shallow water, the 
QSS approximation leads to some improvement over simple, single-scattering 
theories, but it fails to fully account for the effects of multiple scattering. And, 
as used in this paper, it totally neglects multiple reflections between the bottom 
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and the water surface. Novortholoss, tho QSS model provides a framework 
for Identifying potential sources of error In pnBslvo remote measurements of 
bottom depth. 

If a sufficiently homogeneous body of water could be found, it might be 
possible to measure the surface parameters needed for the calculations at one 
(or several) location in the viewed area and to use the remote data to extrapolate 
to other locations. But waters of sufficient homogeneity are relatively rare 
in the world oceans. Even when homogeneous waters can be found, there still 
exists the potential for undetected changes in water turbidity, and hence some 
uncertainty in the depth calculations. 

A conceivable, though not very likely, solution to the problem might be 
the existence of fundamental (though empirical) relationships of biological- 
chemical-optical origin, which can be used to connect the degree of water 
turbidity (or bottom albedo) with observed spectral changes. If such relationships 
could be discovered, they could be used to detect changes in bottom albedo or 
water turbidity remotely without the need for surface detection of these changes. 

Although the existence of these relationships appears to be somewhat 
unlikely, owing to the complexity of the biological, chemical, and optical 
processes involved, nevertheless, an experimental program aimed at thier 
discovery would be a worthwhile task. 
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TABLE I 

MACROSCOPIC MODEL PENETRATION DEPTHS 
CLEAR OCEAN WATER 


Wavelength 

350 

400 

450 

500 

550 

600 

700 

NM 

mmmm 

8.6 

10.8 

12.2 

11.8 

9.2 

3.9 

m 

Meters 

S.S. Penetration Depth 0 

9 

12 

15 

14 

11 

4 

1.8 

Meters 

1 ■ c * 

21.4 

25.7 

26 

21.2 

13.4 

4.4 

1.6 

Meters 

Q.S.S. Penetration Depth 0 i 

13 

18 

22 

21 

14 

6 

1.8 

Meters 

Approximate horizontal 
visibility distance (4 : C) 

34 

1 

43 

49 

47 

37 

16 

6 

Meters 
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Figure- Captions 

Figure 1, A comparison of SS and QSS Models In predictions of tho radiance 
reflectance as a function of layer depth for clear ocoan water at 4G0 nm 
with a zenith sun. 

Figure 2. QSS Model predictions of the radiance reflectance as a function of 
layer depth for clear ocean water at several wavelengths and for a zenith 
sun. 

Figure 3. SS and QSS model penetration depths for cloar ocean water as a 
function of wavelength for a zenith sun, no atmosphere, and a smooth sea 
surface. 

Figure 4. SS and QSS Model radiance reflectance versus layer deptli averaged 
over the two visible-wavelength ERTS-1 MSS bands for clear ocean water 
and a solar zenith angle of 30 degrees. 

Figure 5. SS and QSS Model radiance reflectances versus layer depth averaged 
over the two-visible-wavelength ERTS-1 MSS bands for turbid water and a 
solar zenith angle of 30 degrees. 

Figure 6. QSS Model radiance reflectances versus layer depth averaged over 
the ERTS-1 MSS bands for clear water and several different bottom 
albedos, r Q . 

Figure 7. QSS model radiance reflectances versus layer depth averaged over 

the ERTS-1 MSS bands for turbid water and several different bottom 

albedos, r . 
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Figure 1. A Comparison of SS and QSS Models in Calculations of the Radiance Reflectance as 
a Function of Layer Depth for Clear Ocean Water at 450nm with a Zenith Sun. 
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CLEAR WATER CASE 
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b— W avelength ERTS-1 MSS Bands for Turbid Water 
and a Solar Zenith Angle of 30 Degrees. 
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DEPTH IN METERS 

Figure 6. QSS Model Radiance Reflectances Versus Layer Depth Averaged 
Over the EKTS-1 MSS Bands for Clear Water and 
Several Different Bottom Albedos r 0 . 
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Figure 7. QSS Model Radiance Reflectances Versus Layer Depth 
Averaged Over the ERTS-1 MSS Bands for Turbid Water 
and Several Different Bottom Albedos r 0 . 
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